Magnetic-field-induced binding of few-electron systems in shallow quantum dots 
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Binding of few-electron systems in two-dimensional potential cavities in the presence of an external 
magnetic field is studied with the exact diagonalization approach. We demonstrate that for shallow 
cavities the few-electron system becomes bound only under the application of a strong magnetic 
field. The critical value of the depth of the cavity allowing the formation of a bound state decreases 
with magnetic field in a non-smooth fashion, due to the increasing angular momentum of the first 
bound state. In the high magnetic field limit the binding energies and the critical values for the 
depth of the potential cavity allowing the formation of a bound system tend to the classical values. 

PACS numbers: 73.21. La, 73.43.-f, 85.35.Be 



Quantum dots [jj] (QDs) formed by weak confinement 
potentials are susceptible to perturbations by external 
charge defects and interface roughness 0. The spectra 
of disordered QDs 0, 0| exhibit complex behavior in an 
external magnetic field including pairing and bunching of 
the charging lines, indicating the formation of multiple 
charged puddles in local potential minima inside the large 
QD area. In the quantum Hall effect disorder is held re- 
sponsible for the formation of localized electron reservoirs 
that allows for the pinning of the Fermi level between the 
Landau levels 0. Coulomb blockade of charge confined 
within these tiny reservoirs in quantum Hall systems have 
recently been demonstrated p through scanning probe 
measurements. In this paper we consider electrons local- 
ized in a nanometer-size reservoir created by a potential 
fluctuation. The purpose of the present work is to de- 
termine the role of the magnetic field in the formation of 
bound few-electron states. 

The present study is related to the problem of the ex- 
istence of bound excited states of the negative hydrogen 
ion or charged donor center in a semiconduc- 
tor. In the absence of a magnetic field the H~ has only 
one bound state 0, but an infinite number of weakly 
bound and arbitrarily extended excited states appears 
in weak magnetic field (B — *■ 0+) due to the interplay of 
the magnetic field induced localization and the attractive 
force exerted on the extra electron by the neutral atom 
H, or donor D° This force results from the elec- 
trostatic potential of D'^ which is attractive everywhere 
In two-dimensional quantum wells the D° potential 
for the extra electron is repulsive at large distances 
which prevents the formation of extended excited bound 
states in the weak magnetic field limit. For off-center 
Z)~ in quantum wells magnetic field induces angular mo- 
mentum transitions are found [inj . well-known from the 
theory of QDs States with angular momenta L > 4 
cease to be bound which results in magnetic field induced 



D evaporation [lOj. 

In the symmetric gauge the Hamiltonian of the elec- 
tron in the magnetic field (0, 0, B) contains a diamag- 
netic term Vd{x, y) = muj'^{x'^ + y'^)/8 similar to the two- 
dimensional harmonic oscillator potential with the cy- 
clotron frequency ujc = eBjm. Due to the parabolic form 
of Vd one often assumes [13 that a strong magnetic field 
overcomes the Coulomb repulsion between the electrons 
and makes them orbit around the center of mass even in 
the absence of an external confinement potential. Below 
we demonstrate, that although indeed the magnetic field 
facilitates the binding of few-electron states, even in the 
_B — > 00 limit a finite external potential is necessary to 
stabilize the system. Moreover, we show that the po- 
tential parameters stabilizing the system in the B ^ 00 
limit can be determined using classical physics. 

For the purpose of the present study we need a con- 
finement model realistic enough to account for the un- 
bound states. We choose the QD model of a shallow 
two-dimensional Gaussian well. The Gaussian potential 
is not only the simplest choice but is also realized in verti- 
cal QDs ^3| with small 0| gates. The discussed physics 
is not limited to a Gaussian potential and is independent 
of the specific choice of the confinement provided that it 
has a finite size and depth. A weak magnetic field in- 
duces an infinite number of single-electron states weakly 
bound in the cavity. On the other hand the binding of 
extended two-electron states at weak magnetic fields is 
excluded since the second electron outside the cavity per- 
ceives only the Coulomb repulsion of the electron bound 
in it. In this paper we look for the effect of the magnetic 
field that overcomes the Coulomb blockade and allows 
binding of additional electrons in the cavity. This effect 
can, for example, be used to entangle spins of the 
electrons localized in the cavity and the unbound elec- 
trons. 

Hamiltonian of a single electron with the magnetic 
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field (0, 0, B) perpendicular to the plane of confinement 
reads (symmetric gauge), H — — ?i^V^/2m + Vd{x,y) — 
^ihcOc{xd/dy — yd/dx) + V{x, y), where m stands for the 
effective band mass (GaAs material parameters are used) 
and V{x,y) = — Vbexp(— (x^ +y'^)/R^) is the potential 
well of depth Vq and radius R (for the numerical cal- 
culations we take i? = 50 nm). In the absence of the 
confinement {V = 0) the diamagnetic term Vd{x,y) in- 
duces localization of the single-electron wave functions, 
but by itself it can not keep several electrons together, 
since the Hamiltonian eigenstates can be localized any- 
where. In particular the lowest Landau level (LLL) eigen- 
function (with eigenenergy E^ll = fi-^cl'^) can be put in 
^ (a/27r)i/2exp(-(a/4)((a;-X)2-H 
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the form [If 

{y - Yf) + ieB/2h{xY - yX)), where {X,Y) is the ar- 
bitrary center of the Landau orbit and a = eB /Ti. 
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FIG. 1: (color online) Binding energy of two-electron system 
vs potential well depth for B = Q and 5T. Solid lines were 
calculated using the multicenter basis (the inset shows the 
choice of centers, d is a variational parameter) and the dashed 
lines with the finite-difference scheme. Arrow indicates where 
the ground state angular momentum L changes. 

We studied the dissociation of the two-electron con- 
fined system using the multicenter configuration inter- 
action approach |l7l |. tailored in a way to account for 
both the bound and unbound states of the two-electron 
system. The single-electron states were obtained in the 
basis of the ^'^^'x^' wave functions ^3 in which the po- 
sitions of the centers (X, Y) as well as exponents a are 
treated as variational parameters. The basis generates 
also the excited Landau levels and allows for a precise 
solution of the few-electron Schrodinger equation [l7j . 
We put 15 centers inside the dot (see: inset to Fig. 1) 
and a single center localized far away from the dot (at 
a distance of Z = 2 cm). The method ejects an electron 
from the dot into the LLL localized around the distant 
center when the two-electron dot-confined system is not 
bound. For the unbound N electron states the binding 
energy Eb{N) = En — {En-i + Elll) is not negative. 
Note, that the binding energy is the chemical potential 
fiN = En — En -I of the iV-electron system determining 



the single-electron charging of the QD devices calculated 
with respect to the LLL Eb{B) = hn — Elll- 

The calculated binding energies are plotted in Fig. 1 
by the thick solid lines as function of the depth of the 
well. For B = the two-electron system becomes bound 
in the L = state near Vq = 5.75 meV. For 5T binding 
appears for a shallower well with Vq = 4.2 meV, and the 
first bound state has L = 4. 

To verify these results, we used an alternative approach 
|l8j using single-electron wave functions diagonalized nu- 
merically on a finite-difference mesh. The method [l^ 
allows for an exact description of states localized in a cir- 
cular potential of arbitrary profile, but does not account 
for the delocalized LLL's since the calculations are per- 
formed within a box of finite size. The results are shown 
by the dashed lines in Fig. 1 and agree very well with 
our previous approach. Above the dissociation threshold 
Eb ~ the dashed lines become positive due to the ab- 
sence of the distant LLL in the basis The critical 
Vq value is well reproduced. For B = T, the finite dif- 
ference approach works better than the multicenter basis 
(see Fig. 1), because it better describes the tails of the 
wave functions for weakly bound states, which in the ab- 
sence of a magnetic field can be arbitrarily long. In the 
following we give the results as obtained with the finite 
difference method [l^ . 
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FIG. 2: (color online) Binding energy of the two-electron sys- 
tem vs the B for depth of the potential from Vo ~ 7 meV 
(lowest curve) to 4.18 meV (highest curve). Horizontal dashed 
lines show the binding energies of the classical system. The 
energy cusps due to the angular momentum transition are 
linked by the thin solid lines. The dotted curve shows the 
binding energy of the state with L — 3 for Vo = 5.4 meV. 

The two-electron binding energy as function of the 
magnetic field is plotted in Fig. 2 for various depths 
of the dot. The two electrons are bound for any B when 
Vb = 6 meV, or more. The bound ground state under- 
goes interaction-induced angular momentum transitions 
|ll| . The interaction energy increases with increasing 
magnetic field till it reaches a value at which the transi- 
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tion to a more weakly localized state of lower interaction 
energy is more favorable. At high magnetic fields the 
angular momentum transitions keep the size of the 
system close to the size of its classical counterpart [l^ . 
The angular momentum transitions are associated with 
the ground-state spin oscillations. States with even L 
correspond to spin singlets (5 = 0), and those with odd 
L to spin triplets {S = 1). 

For Vo = 5.4 meV bound states appear only in the 
presence of an external magnetic field (see Fig. 2). Notice 
that, for shallower QDs the first bound state has non- 
zero L. Wc found that states of specific L arc bound 
only within a finite interval of the magnetic field (see the 
binding energy for i = 3 at Vq = 5.4 mcV plotted by the 
dashed line). 
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FIG. 3: (color online) Binding energy of the two-electron sys- 
tem minus the classical binding energy vs the inverse of the 
magnetic field for depth of the potential from Vo ~ 7 meV 
(highest curve) to 4.18 meV (lowest curve). The thin dashed 
straight lines are given by 0.39/B , 0.78/B and 1.28/B. 

The horizonal dashed lines in Fig. 2 show the val- 
ues of the binding energies obtained for classical point 
charges |l9l |2C| , via minimization of the two-electron po- 
tential energy ~ j^^^~2V{'w/2,0) with respect 
to the interclcctron distance w. The quantum binding 
energies tend to the classical values at high B. This is 
shown clearly when wc plot the difference between the 
quantum and classical values as function of the inverse 
of the magnetic field (sec Fig. 3). Wc conclude that the 
binding energy tend as 1 /B to the classical value in the 
limit of high magnetic fields The classical charac- 
ter of the binding is due to the point localization of the 
single-electron LLL wave functions in the high magnetic 
field. 

The cusps of the three-electron binding energy are due 
to both the TV = 3 and N ~ 2 angular momentum tran- 
sitions (cf. Fig. 4) . The 1/B convergence of the binding 
energy to the classical limit is also found for three elec- 
trons (see Fig. 5). 

Fig. 6 shows the 'phase diagram' for the maximum 
number of bound electrons as function of the well depth 




FIG. 4: (color online) Binding energy of the three-electron 
system vs the magnetic field for three values of the poten- 
tial depth. Horizontal dashed lines show the binding energies 
of the classical system. The energy cusps corresponding to 
the angular momentum transitions in the three- (two-) elec- 
tron ground state are linked by thin solid (dotted) lines. The 
numbers below the lowest curve list L for N = 3. 
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FIG. 5: (color online) Same as Fig. 3 but for = 3. The thin 
dashed straight lines are given by 6.6/B, 5.4/_B and 4.1/_B. 



and the magnetic field. The thin lines joining the cusps 
of the phase diagram boundaries corresponding to iV = 2 
angular momentum transitions. The cusps related to the 
iV = 3 angular-momentum transitions are marked with 
arrows. Dashed lines show the depth of the potential 
well which allows binding of the second and third classi- 
cal electron . In Fig. 7 we plot the phase diagram as 
function of \/B to demonstrate the convergence of the 
phase boundaries to the classical result at infinite mag- 
netic field. 

In an external magnetic field the diamagnetic effect in- 
creases faster the free-electron LLL than the energies of 
the states localized in the cavity. This leads to the energy 
shift of the ionization threshold observed in a magnetic 
field. According to our results the magnetic field induces 
the formation of bound electron states in QDs. Such a 
formation was actually observed in single-electron charg- 
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FIG. 6: (color online) The maximum number of bound elec- 
trons on the Vo/B plane. Below the lowest curve only one 
electron can be bound, and above the highest the third elec- 
tron becomes bound, between the curves the maximum num- 
ber of bound electrons is two. The horizontal dashed lines 
show the classical depths of the well allowing binding of the 
second and third electron. 
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FIG. 7: (color online) Same as Fig. 6 but in 1/B scale. The 
thin dashed straight lines are given by V^'"'"" + An/B, with 

yclass ^ 4 j^gy^ yclass ^ g 345 ^ ^ Q.T, As = 5. 

ing experiments on self-assembled QDs |23|. The QDs 
[231 at low magnetic fields bind six electrons. Only for 
B > 8 T the capacitance signal related to the binding of 
sixth and seventh electron emerge from the background 
of the LLL wetting layer charging (see Fig. 3 of Ref.|2^). 
Self-assembled QDs have a fixed depth of the confinement 
potential, but the phase diagram (Figs. 6 and 7) can be 
experimentally determined in the capacitance measure- 
ment of the vertical quantum dot device . The shape 
of the confinement potential is exactly Gaussian for a 
proper choice of the size of the gate (see Fig. 6 of Q) 
forming a small electrostatic QD. The depth of the poten- 
tial is simply changed by the gate voltage. The number of 
bound electrons can be determined for any magnetic field 
by counting the capacitance peaks appearing as function 
of the gate voltage below the two-dimensional electron 
gas charging (like the one of Fig. 3 Ref. H^l). 

The discussed effect can be used to entangle electron 
spins [isl . One needs a cavity which at i? = binds only 
a single electron. The second electron is bound by the 



field and one can arrange for the bound state to have zero 
spin (even L). After the magnetic field is switched off the 
second electron will become unbound and is repelled by 
the dot and thus goes with its spin state entangled with 
the localized electron. 

In summary, we have shown that a strong magnetic 
field induces the formation of bound few-electron states 
in potential cavities that are deep enough to bind the 
classical point electrons. 
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